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ABSTRACT 

We present results from an XMM-Newton observation of the nearby Seyfert 2 galaxy 
NGC 5643. The nucleus exhibits a very flat X-ray continuum above 2 keV, together 
with a prominent K^, fluorescent iron line. This indicates heavy obscuration. We mea- 
sure an absorbing column density Nh in the range 6-10x10^^ cm^^, cither directly 
covering the nuclear emission, or covering its Compton-reflection. In the latter case, 
we might be observing a rather unusual geometry for the absorber, whereby reflection 
from the inner far side of a torus is in turn obscured by its near side outer atmosphere. 
The nuclear emission might be then either covered by a Compton-thick absorber, or 
undergoing a transient state of low activity. A second source (christened "X-1" in this 
paper) at the outskirts of NGC 5643 optical surface outshines the nucleus in X-rays. If 
belonging to NGC 5643, it is the third brightest {Lx ~ 4 x 10^" erg s~^) known Ultra 
Luminous X-ray source. Comparison with past large aperture spectra of NGC 5643 
unveils dramatic X-ray spectral changes above 1 keV. We interpret them as due to 
variability of the active nucleus and of source X-1 intrinsic X-ray powers by a factor 
>10 and 5, respectively. 

Key words: galaxies: active - galaxies individual (NGC 5643) - galaxies muclei - 
galaxies: Seyfert - X-rays:galaxies 



1 INTRODUCTION 

The X-ray spectra of Seyfert 2 galaxies exhibit significant 
photoelectric absorption (Warwick et al. 1989; Awaki et al. 
1991; Turner et al. 1997; Risaliti et al. 1999). This evi- 
dence has been interpreted as supporting the predictions 
of the Seyfert unification scenarios, whereby "type 2" ob- 
jects are seen at high inclination angles with respect to an 
azimuthally-symmetric gas and dust structure (the "torus" ) , 
which prevents us from directly observing the nucleus and 
the Broad Line Regions. ASCA observations, however, un- 
veiled Seyfert 2s with peculiar X-ray spectral properties: 
an inverted spectrum (energy index, a < 0); large Equiva- 
lent Width (EW) Ka iron fluorescent lines (EW from a few 
hundreds to thousands eV); very little or no evidence for 
photoelectric absorption (Elvis & Lawrence 1988; Koyama 
et al. 1989; Ueno et al. 1994). This phenomenology was 
interpreted as due to the column density of the absorber 
in this subclass of Seyfert 2s being so large, that it sup- 
presses almost to invisibility the X-ray continuum, and en- 
hances the contrast between the iron feature and the un- 
derlying continuum. When the absorber column density ex- 



ceeds Nh — cTj"^ ~ 1-5 X 10^* cm~'^, the absorbing matter 
is optically thick to Compton scattering, and the nuclear 
photons are downscattered to energies where the photoelec- 
tric absorbing cross section becomes dominant. The nuclear 
continuum is then substantially suppressed below 10 keV, 
and the nucleus can be seen only along reflected and/or 
scattered optical paths, which do not intercept the obscur- 
ing matter. Indeed, later BeppoSAX observations, which ex- 
tended the sensitive bandpass beyond 10 keV, confirmed this 
hypothesis by detecting nuclear emission piercing through 
<; 2 X 10^* cm"^ absorbers (Matt et al. 1997b; Vignati et 
al. 1999; Guainazzi et al. 2000). All the above justify the 
nomenclature of "Compton-thick" given to this sub-class of 
Seyfert 2s. 

The X-ray spectra of the three closest AGN (Centau- 
rus A, NGC 4945, and the Circinus Galaxy) are absorbed 
by column densities larger than 10^^ cm~^, and two of them 
are Compton-thick. This simple fact suggests that Compton- 
thick Seyfert 2 galaxies may constitute a large fraction of 
their parent population. Indeed, BeppoSAX measurements 
of large samples of optically-selected Seyfert 2 galaxies sug- 
gested that they may constitute a fraction as large as 50 per 
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Table 1. Log of the observations discTisscd in tliis paper 



Observatory-instrument /s 

ASCA-SIS/GIS 
BeppoSAX-LECS/MECS 
ROSAT-HRI 
XMM-Newton-MOS/pn 



Date 

21-Feb-1996 
Ol-Mar-1997 
28- Aug- 1997 
08-Feb-2002 



Exposure time 
(ks) 
36.7/41.5 
7.1/10.4 

10.0 
9.4/7.1 



the source reference frame; errors arc at the 90 per cent 
confidence level for 1 interesting parameter for the best-fit 
model parameters, and at the 1-a level for count rates; chem- 
ical mixture follows the Anders & Grevesse (1989) measure- 
ment; H„ = 70 km Mpc"^ s"^ (Bennett et al. 2003). At the 
NGC 5643 distance (16.9 Mpc), 1" corresponds to 82 pc. 



cent of the total population (Maiolino et al. 1998, Risaliti 
et al. 1999), at least in the local universe. It is hard to tell 
whether the fraction of heavily obscured objects remains 
that large at cosmological distances. If so, this may have 
important consequences for the history of accretion, and for 
the total energy budget ratio between accretion and stellar 
light in the universe (Fabian 1999). 

We have started a program to study a complete, unbi- 
ased, optically-defined sample of Compton-thick Seyfert 2 
galaxies (Risaliti et al. 1999) with XMM-Newton (Jansen 
et al. 2001), with the main goal of characterizing their X- 
ray spectral properties in the hard X-ray regime. Prelimi- 
nary results on the whole sample arc discussed by Guainazzi 
et al. (2004). In this paper, we present the observation of 
NGC 5643. It allows us to describe in details the analysis 
method followed in the study of the sample. The interest of 
this specific observations is two-fold: a) the large collecting 
area of the XMM-Ncwton optics allows us to measure a, al- 
beit extreme, Compton-thin absorber; b) better imaging res- 
olution with respect to previous mission allows us to discover 
a serendipitous source, which apparently belongs to the op- 
tical/UV surface of the host galaxy, and outshines the nu- 
cleus. If this source is indeed associated with the NGC 5643 
host galaxy, it represent one of the brightest Ultra-Luminous 
X-ray (ULX) source ever observed, with a total X-ray lumi- 
nosity larger then 4 x 10*° erg s~^. 

1.1 Some properties of NGC 5643 

NGC 5643 is a nearby {z = 0.004) SAB(rs)C galaxy, known 
to host a low- luminosity Scyfcrt 2 nucleus (Phillips et al. 
1983). It exhibits an extended emission line region elongated 
in a direction close to the radio position angle (Morris et 
al. 1985), due to a "V" -shaped structure of highly excited 
gas. This is probably the projection of a 1.8 kpc, one-sided 
ionization cone (Simpson et al. 1997). Perpendicular to the 
radio axis, a dust lane covers the nucleus (Simpson et al. 

1997) . The fact that NGC 5643 belongs to the class of "ex- 
treme infrared" galaxies (Antonucci & Olszewski 1985) at- 
tracted rrmch attention in the past to try to explain the ori- 
gin of IR emission in this Active Galactic Nucleus (AGN). 
Although intense episodes of star formation are occurring 
in its nearly circular arms, evidence with respect to the nu- 
cleus is still controversial. Mid-IR diagnostics suggest that 
the AGN dominates the IR energy budget (Genzel et al. 

1998) . Comparison of optical spectra with synthesis models 
are, however, consistent with a "starburst/Seyfert 2 compos- 
ite" spectrum (Cid Fernandes et al. 2001). The NGC 5643 
nucleus is a strong radio emitter as well, most likely powered 
by the AGN (Kewley et al. 2000). 

The log of the observations discussed is this paper is 
reported in Table 1. In this paper: energies are quoted in 



2 XMM-NEWTON OBSERVATION 

XMM-Newton (Jansen et al. 2001) observed the sky re- 
gion encompassing NGC 5643 on February 8, 2003. The X- 
ray imaging EPIC cameras (pn, Striider et al. 2001; MOS, 
Turner et al. 2001) were operated in Full Frame Mode, with 
the Medium optical photons blocking filter. The Optical 
Monitor (OM; Mason et al. 2001) was operated in standard 
Image Mode with the UVWl filter, sensitive in the 2500- 
4000A bandpass. Data reduction was performed with ver- 
sion 5.4.1 of the SAS software package (Jansen et al. 2001). 
Wo employed the most updated calibration files available at 
the time the reduction was performed (June 2003). None of 
the X-ray sources discussed in this paper was detected by 
the high-resolution spectroscopy camera, RGS (der Herder 
et al. 2001). Standard data screening criteria were applied 
in the extraction of scientific products. The particle back- 
ground stayed at a quiescent level (Lumb et al. 2002) during 
the whole observation, making unnecessary any background- 
rejection filtering. Pattern to 4 (12) were employed in the 
extraction of pn (MOS) scientific products. Spectra and light 
curves were extracted from circular regions of 25" radius 
around the source centroids. Response matrices, appropriate 
for each of the spectra discussed in this paper, were gener- 
ated with the SAS tasks arfgen and rmfgen. Background 
spectra were extracted from source-free regions of the same 
CCD as where the source is located. Spectra were rebinned 
to ensure that: a) the instrumental energy resolution is over- 
sampled by a factor not larger than 3; b) at least 25 (50) 
counts are present in each MOS (pn) spectral channel, to 
ensure the applicability of the statistics in the evaluation 
of the fit quality. Images, spectra and response matrices of 
the two MOS cameras were combined together, to increase 
signal-to-noise. Spectral fits were performed simultaneously 
on both EPIC cameras, using the 0.35 12 and 0.5 10 kcV 
energy bands for the pn and the MOS, respectively. As no 
significant variability was detected during the XMM-Newton 
observation, in this paper we will discuss the time-averaged 
spectra only. 

In Fig. 1 the 0.5-10 keV sky coordinates MOS image 
of the innermost 3' around the NGC 5643 nucleus is shown. 
Four sources are detected by the EPIC cameras in this field 
at a signal-to-noise ratio larger than 3. The source detected 
at (Q2000 = 14'*32'"40.9", ^2000 = -44°10'26") coincides 
with NGC 5643 optical nucleus within the typical XMM- 
Newton attitude reconstruction accuracy. Its total counts 
are 1380 ± 40 and 640 ± 30 in the pn and MOS camera, 
respectively. The source located ~0.8'NE the rmclcus is 50 
per cent brighter then the nucleus. We will refer to it as 
"NGC 5643 X-1" hereafter. Its total counts are 2030 ± 60 
and 1320 ± 40 in the pn and MOS cameras, respectively. 
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NGC5643 
XMM-Newton MOS 




-44°09'00' 



Figure 1. 0.5—10 keV single- and double-events MOS image, re- 
stricted to the 3' around the NGC 5643 optical nucleus. The image 
is smoothed with a Gaussian function (cr = 1.5"). The scale in is 
units of unsmoothcd pixel counts. 




Energy (keV) 

mguainaz 9-Feb-20i 

Figure 2. NGC 5643 nucleus EPIC spectra (upper panel; dots: 
pn; crosses: MOS) and residuals in units of standard deviations 
(lower panel) against a power-law model. In the inset: data/model 
ratio in the 5.5-7.5 keV energy band against a power-law model 
above 2 keV (pn: filled circles; MOS empty squares). In the latter 
panel, data are rebinned according to a constant factor SE = 
50 eV for plotting purposes only 



2.1 Spectral analysis of the nucleus 

In Fig. 2 we show the resuh of a power-law fit* on the 
NGC 5643 nucleus spectrum. The fit is clearly unacceptable 
(x^/*^ = 417.0/29). As often observed in highly absorbed 



Seyfert 2 galaxies (Matt et al. 2000) the residuals exhibit a 
"turning point" around 2 keV, suggesting that the contin- 
uum needs to be modeled with (at least) two components. 
Additionally, a broad excess emission feature is localized 
around E ~ 0.85 keV (observer's frame), together with a 
narrow emission feature around 6 keV. In order to achieve 
a full physical characterization of the overall spectrum, we 
have first separately analyzed the regimes above and below 
the spectral "turning point" . 



2.1.1 The soft X-ray spectrum (E < 2 keV) 

The soft energy spectrum is remarkably complex. A fit with 
a single featureless continuum in the energy band redwards 
of 1 keV is clearly inadequate (e.g.: xt — 3.5 if a power-law is 
employed), and leaves large positive residuals in the energy 
range between 0.8 and 1 keV. In this region a "forest" of 
iron-L transitions exists. At least two thermal components^^ 
are required in order to adequately account for this feature, 
with temperatures kTs ~ 0.17 keV and kTh ~ 0.69 keV 
ix'^/'^ = 9.7/8). Alternative "two continuum components" 
descriptions of the soft X-ray spectrum are possible, and 
yield statistically equivalent fit, such as the combination of 
thermal emission and power-law (kT ~ 1.1 keV; F ~ 2.2), 
or of thermal and multi-temperature accretion disk black- 
body emission (model diskbb in XspeC; kT*'"'''""'' ~ 90 eV; 
kT'^'^'' ~ 1.0 keV). 

Is it alternatively possible that the soft X-ray spectrum 
is dominated by scattering? The soft X-ray spectrum can 
be in principle reasonably fit with a combination of a, say, 
power-law continuum plus three emission lines (xt — 1.23). 
However the spectral index in this scenario (F ~ 4.0) is 
too steep to represent the mirror of the AGN intrinsic 
continuum, even if self absorption effect in the scattering 
plasma are neglected. However, high-resolution grating mea- 
surements of bright absorbed Seyfert galaxies are showing 
that the soft X-ray emission is often dominated by heavily 
blended emission lines (Kinkhabwala et al. 2002; Brinkmann 
et al. 2002; Sako et al 2000). The blend can mimic continuum 
emission - alongside a few emission "peaks" corresponding 
to the brightest well defined lines - in low-resolution spec- 
tra (Guainazzi et al. 1999; Cappi et al. 1999). Testing this 
hypothesis would require high-resolution good quality spec- 
troscopic data, which are not available for the NGC 5643 
nucleus. Nonetheless, we fit the EPIC spectra with a com- 
bination of emission lines only, aiming at inferring the main 
qualitative properties of the line emitting plasma in this 
scenario. In order to avoid line blurring effects induced by 
spectral rebinning, we fit the unbinned spectrum, using the 
C-statistics (Cash 1976). Only the pn spectrum was used, 
due to larger effective area. No background subtraction was 
performed, as its contribution is < 10 per cent in the soft 
X-ray band. Seven lines are required to fit the spectrum (see 
Fig. 3, and Table 2). Typical statistical uncertainties on their 
centroid energy are ~ 10 eV (except for the highest energy 



* All the models employed in this paper include photoelectric 
absorption by a column density due to the contribution of neu- 
tral matter in our Galaxy along the line-of-sight to NGC 5643: 
NH,Gal = 8.3 X 10^0 cm-2 (Dickey & Lockman 1990) 



t hereafter thermal emission from an optically thin, coUisionally 
excited plasma is modeled via the code mekal in XsPEC (Mewe et 
al. 1985). Unless otherwise specified, the plasma abundance has 
been kept fixed to Z = 0.5Zq 
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X (Angstrom) 



Figure 3. 10-35A pn spectrum (crosses) superposed to the 
emission-line only model best-fit {solid line - details in text) 



Table 2. Centroid energies for the measured emission lines, if 
the soft X-ray pn spectrum of the NGC 5643 nucleus is fit with a 
model constituted by emission lines only. The right column shown 
the possible identifications. They take into account the statisti- 
cal uncertainties on the centroid energies only. In brackets, ad- 
ditional identifications are listed when systcmatics uncertainties 
are taken into account as well. RRC = Radiative Recombination 
Continuum 



Energy (keV) Identifications 

0.428 Cvi /? (Nvi) 

0.566 Ovii 

0.668 Ovii 7 

0.787 Fcxvii 

0.899 Neix, Fexix 

1.028 Neix 13, Fexxii 

1.176 Neix RRC, Fexxiv 



line in Table 2, for which they are ~ 20 eV). Systematic 
uncertainties on the gain are of the same order. 

Better data quality is necessary to achieve a fully un- 
ambiguous dcconvolutiori of the soft X-ray spectrum. For 
sake of simplicity, we will assume in the following a model- 
ing of the soft X-ray spectrum in terms of a two-component 
optically-thin thermal plasma. 

2.1.2 The hard X-ray spectrum (E > 2 keV) 

The nuclear spectrum above 2 keV is no less complex. A 
formally adequate fit is achieved through an unabsorbed 
power-law and an emission line (parameterized through a 
simple Gaussian profile). The latter is required at the 99.87 
per cent confidence level according to the F-test. However, 
the best-fit parameters suggest that this is no more than 
a purely mathematical description. The power-law spectral 
index is inverted (photon index, F ~ —0.8), what is totally 
unusual for AGN (Nandra ct al. 1997; et al. 1997; Reeves & 
Turner 2000; Perola et al. 2002). The emission line centroid 
energy {Ec ~ 6.42 keV) is consistent with fluorescent 
emission from neutral or mildly ionized iron. The formal 



line Equivalent Width is huge: EW ~ 2.8 keV (cf. the inset 
in Fig. 2). 

Flat hard X-ray spectra with huge EW fluorescent 
iron lines in Seyfert 2 galaxies can be produced if the bulk 
of hard X-ray emission is due to Compton-reflection of an 
otherwise invisible "standard" AGN nuclear continuum, ei- 
ther because the AGN switched to an "off" state during the 
XMM-Ncwton observation, or because its emission is totally 
suppressed by a Compton-thick absorber (Matt et al. 2003). 
We have therefore fit the EPIC nuclear spectra in their full 
energy bandpasses with a combination of two thermal com- 
ponents, a Gaussian emission line, and a "bare" Compton- 
reflection component (model pexrav in XSPEC; Magdziarz & 
Zdziarski 1995). The fit is unacceptable (xV^ = 49.9/22), 
mainly due to large excess residuals between 2 and 5 keV. 
This excess could be in principle due to optically thin elec- 
tron scattering of the otherwise invisible nuclear continuum. 
Indeed, adding a power-law component yields a statistically 
acceptable fit {x^ /v = 18.6/21). However, the spectral index 
of this power-law is too flat (F = — O.lioie)! for its interpre- 
tation as the "warm mirror" of the nuclear continuum to 
be plausible. The combination of a Compton-thick absorber 
with nuclear electron scattering can be therefore ruled out. 

A possible alternative is that the Compton-reflection 
component is in turn absorbed by matter with a covering 
fraction lower than 1. In this case, the 2-5 keV spectrum 
is accounted for by ~5% of the Compton-reflection contin- 
uum (see Model#l in Table 3) which "leaks" through its 
absorber. This leakage has the appropriate flatness. An ab- 
sorber with a ~95 per cent covering fraction and Nh = 
(1.0 ± 0.3) X 10'^'* cm~^ is required in this scenario. If the 
hard X-ray continuum is dominated by Compton-reflection, 
Kq, emission lines from O to Cr are expected to imprint their 
signatures along with the Fe line. They are not detected in 
the EPIC spectra of NGC 5643. However, their EW upper 
limits (200-300) eV are inconclusive (Matt et al. 1997b). 

The depth of the K photoelectric absorption edge is 
in principle a powerful diagnostic for the metallicity of 
the absorbing and/or reflecting matter. A further 97 per 
cent confldence level improvement in the quality of the fit 
(Ax^/Av = 5.0/1) is obtained in Model#l of Table 3 if 
the iron metallicity is left free. If this change is attributed 
to one of the two components only (what is unlikely, if ab- 
sorber and reflector represent one and the same system; see 
Sect. 4.1): Zpe < Q.7 Zq. 

Alternatively, Compton-thin obscuration - still allowing 
a fraction of the direct nuclear emission to pierce through 
the X-ray absorber in the EPIC energy bandpass - can 
as well significantly harden a typical AGN spectrum. We 
have therefore substituted the partially covered Compton- 
reflection component with a power-law modified by pho- 
toelectric absorption. Again the fit is statistically accept- 
able {x't — 0.95). The best- fit models results are reported 
in Table 3. In this scenario the 2-5 keV excess can be ac- 
counted for by optically-thin scattering of the nuclear con- 
tinuum, with a reasonable spectral shape (F = 1.4±o:2)- 
The extrapolation of the scattered nuclear continuum into 
the 0.4-1.2 keV band is ~15 per cent in Model #3. No sig- 
nificant evidence is found for a deviation of the absorber 
metallicity from standard solar abundances in this scenario: 
Zpe = 0.8±o.2 Ze- In Model#2 of Table 3 the electron scat- 
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Table 3. Best-fit parameters and results for the EPIC NGC 5643 
nuclear spectrum. Model #1 is constituted by the combination of 
two thermal components (with temperatures kTs, and feT/j), a 
Compton-reflection continuum, partly absorbed by matter with 

a covering fraction fc and column density Nu, and a Gaussian 
emission line profile with centroid energy Ec, intensity lune, and 
Equivalent Width EW . Model #2 is constituted by the combina- 
tion of three thermal components (the hottest with temperature 
kT^f^), an absorbed power-law, and the Gaussian emission line 
profile. In Model #3 the "ultra-hard" thermal component is sub- 
stituted by a locally unabsorbed power-law. fs is defined as the 
scattering fraction between the transmitted and the scattered nu- 
clear component in Model #3 



Parameter 


Model #1 


Model #2 


Model #3 


r 


2.3±1i 


1 6±i° 




Nh (10^3 cm-2) 


10 ±3 


6±4 


7±2 


fc or fs (per cent) 


94.7±3-| 






Ec (kcV) 


6-39±0:g| 


c 4Q+o.a2 


a 4Q+o:02 


r a 

-'-line 


8±4 


1.4 ±0.4 


i5±8-2' 

^•'^='=0.3 


EW (eV) 


760 


490 


500 


kTs (keV) 


n 17+0.07 


n 1^+0.03 


n 1 c+0.03 


kn (keV) 


n 7Q_i_0.10 


0.67 ± 0.05 


n fi7+0.05 
"•"'='=0.04 


kT^h (keV) 




> 5 




0.5-2 keV flux'' 






2.16 ±0.15 


2-10 keV flux*' 






8.4±}:3 




21.2/19 


20.0/21 


20.8/22 



"in units of 10 ® cm ^ s ^ 

''corrected for Galactic absorption, in units of 10"^^ erg cm~^ 



tering is substituted by a "ultra-hot" (fcT ~ 5 keV) thermal 
component. 

In order to investigate the details of the iron line struc- 
ture, we have performed a Cash statistic fit of the unbinned 
pn spectrum in the 4-8 keV energy range, again to prevent 
the line profile from being blurred by the applicability 
requirement. The contribution of the background has been 
estimated by fitting the background spectrum in the same 
energy band with a simple power-law. The best-fit param- 
eters for the background are: T = 3.08 and 1 keV normal- 
ization, N = 5.8 X 10~® cm~^ s~^. The source spectrum 
continuum was than fitted with a double power-law, one 
of the two components having its parameters fixed to the 
background best-fit model. Two lines are detected in the 
spectrum, at a confidence level larger than 99 per cent for 
1 interesting parameter (sec Fig. 4; Lampton et al. 1976). 
They have centroid energies Ei = 6.40 ±0.02 keV and E2 = 
6.93 ±0.05 kcV, and intensities (1.4±g:|) x lO"'^ cm'^ s"\ 
and (4±2) x 10~® cm~^ s~^, respectively. The fainter line is 
consistent with Ka fluorescence of H-like iron. 

The observed fluxes - calculated according to Model 
#3 best-fit parameters in Table 3 - are 1.7 x 10^^^ and 
8.3 X 10"" erg cm"^ s"^ in the 0.5-2 keV and 2-10 keV en- 
ergy ranges, respectively. The corresponding luminosity in 
the 0.5-2 keV energy range, once corrected for the Galac- 
tic absorption, is ~ 9 X 10^® erg s~^. The intrinsic AGN 
luminosity, further corrected for its local absorption, is 
(2.5±^:o) X 10"^ erg s"^ in the 0.5-10 keV energy range. 



2.2 Spectral analysis of NGC 5643 X-1 

In Fig. 5 we show the EPIC spectra of NGC 5643 X-1 and its 




6.5 7 
Energy (keV) 

Figure 4. Upper panel: pn spectrum (crosses) and best-fit model 
(solid line) in the 4—8 keV energy band, showing two emission 
lines with centroid energies ~ 6.40 keV and ~ 6.93 keV. Lower 
panel: residuals against the best-fit model. A constant SE = 50 eV 
rebinning has been applied for plotting purposes only 
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Figure 5. Spectra (upper panel) and residuals in units of stan- 
dard deviation (lower panel), against a photoelectrically absorbed 
power-law for NGC 5643 X-1 



best-fit with a photoelectrically absorbed power-law model. 
This fit is acceptable (^ jv = 46.4/48). An equally good fit 
can be obtained with a thermal brcmsstrahlung continuum. 
A fit with a multitemperature blackbody (diskbb in XSPEC) 
is instead poor. A summary of the best-fit parameters and 
results is shown in Tabic 4. No iron fluorescent emission 
line (or any other narrow band features) is detected in the 
spectrum. The upper limit on the EW of a "narrow" (i.e. 
cr = 0), neutral (E^ = 6.4 keV in the observer's frame) line 
is 270 eV. 

The observed flux in the 0.5-10 keV energy band is 
(1.11 ± 0.06) X 10~^^ erg cm~^ . It corresponds to an 
intrinsic luminosity of (4.4 ±0.2) x lO'*" erg s"^ if NGC 5643 
X-1 is located at the galaxy redshift. 
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Table 4. Best- fit parameters and results for the NGC 5643 X-1 
EPIC spectra. Models legenda: po = power-law; dlskbb = multi- 

temporatnro disk })lack})0(h': brems = thermal })rcmsstrahlung, 



Parameter 
iVjf(102l cm-2; 

r 

kT (keV) 
0.5-2 keV flux" 
2-10 keV flux" 



po 
1.7 ±0.4 

1 fiQ+0.09 



4.4 ±0.5 

7.5 ±0.6 
46.3/48 



dlskbb 
<0.05 



1 42+0.08 



brems 
1 l±0-2 



8.5±2;2 



87.8/48 45.9/48 



° absorption-corrected, in units of 10 ^"^ erg cm ^ s 



Table 5. Sources detected in the ROSAT/HRI observation of the 
NGC 5643 field, within 3' the optical nucleus at a signal-to-noise 
ratio larger than 3 



Table 6. Best-fit parameters and results for the ASCA and Bep- 
poSAX spectra of the 3' radius region encompassing NGC 5643. 
The best-fit continua are: mekal plus power-law for ASCA, and 
power-law for BeppoSAX, respectively. Both models include a 
narrow Gaussian profile, and absorption by a column density 
NH,Gal = 8.3 X lO^u cm-2. 



Parameter 
F 

kT (keV) 
Ec (keV) 
EW (keV) 
0.5-2 keV flux"" 
2-10 keV flux" 



ASCA 
1.04 ±0.14 

0.68±O:i» 
6.45 ± 0.11 

1.7 ±0.4 
o.o=n„ 4 

9.8 ± 1.5 
131.0/184 



BeppoSAX 

1 

^•^^0.3 



6 43±"'i" 

U. lO^o 

5.4±}'| 
10 ±2 
9.1/13 



"absorption-corrected, in units of 10 erg cm ^ s 



a2000 <52()()() Count rate Flux" 
(10--^ s-i) 

14''32"'40.7'' -44°10'24" 7.1 ± 1.0 2.3 ±0.3 

14''32'"41.9^ -44°09'36" 3.2 ± 0.7 1.5 ±0.3 



"0.5-2 keV flux in units of 10~^^ erg cm"^ s~^, absorption- 
corrected and extrapolated from the measured count rate via the 
HEASARC on-line tool Pimms, using the best-fit models of the 
XMM-Newton observation 

3 THE X-RAY HISTORY OF NGC 5643 

In this Section we compare the results of the XMM-Newton 
observation of the NGC 5643 field with earlier ASCA, Bep- 
poSAX, and ROSAT/HRI observations. Data were retrieved 
from the NASA/HEASARC and ASI/ASDC archival facili- 
ties, and reduced according to standard procedures. In the 
HRI image both NGC 5643 XMM-Newton sources are in- 
dividually detected, thanks to its good spatial resolution 
(see Table 5) On the other hand, the broad Point Spread 
Function of the ASCA and BeppoSAX optics mixes irreme- 
diably the contribution of the two sources in their typical 
aperture. In both latter cases, we analyzed a single spec- 
trum extracted from a circular region of 3' radius (except for 
the BeppoSAX LEGS detector, for which an aperture of 2' 
was used due to the lack of appropriate calibrations) . Spec- 
tra were rebinned according to the same criteria employed 
for the XMM-Newton/MOS spectra (cf. Sect. 2.1). Back- 
ground spectra were extracted from nearby regions in the 
same field-of-view where the source is located. The soft flux 
exhibits a dynamical range of about a factor 2 between the 
ASCA (faintest) and the XMM-Newton (brightest) state. 
The variation amplitude in the hard X-rays is somewhat 
smaller (~ 30 per cent). 

Is it possible to unambiguously ascribe the observed 
historical flux variability either to the NGC 5643 nucleus or 
to source X-1? The fainter states caught by BeppoSAX and 
- above all - ASCA may represent a "diluted" version of a 
more dramatic variability afi'ecting only one of the members 
of the pair. Indeed, NGC 5643 X-1 underwent a factor of 
3 variation in flux between the 1997 ROSAT/HRI and the 
2002 XMM-Newton observations. However, a quantitative 
assessment of this issue across the whole X-ray history of 
NGC 5643 is not straightforward. 

The ASCA and BeppoSAX spectra exhibit mutually 



different spectral shapes, notwithstanding the universal 
presence of the large EW Kq, fluorescent iron line (see 
Fig. 6). A simple power- law continuum is an adequate rep- 
resentation of the BeppoSAX spectrum, with a steep index 
(r = 1.9±o;3). This is markedly discrepant with the photon 
index measured by XMM-Newton, when a simple power- 
law is applied above 2 keV (F ~ -0.8; cf. Sect. 2.1.2). The 
ASCA spectra requires at least two continuum components. 
Adopting for simplicity the combination of a thermal and a 
power- law component, the temperature of the former is con- 
sistent with kTh as measured by XMM-Newton, whereas the 
latter's spectral index is intermediate between BeppoSAX 
and XMM-Newton ones. A summary of the ASCA and Bep- 
poSAX best-fits - which represent for each case only one of 
the possible statistically equivalent solutions - is reported in 
Table 6. 

In order to quantitatively estimate the differences with 
respect to the XMM-Newton measurements, we have built 
a "composite" XMM-Newton spectrum, combining the best 
flt models for the NGC 5643 nucleus and source X-1. This 
"composite" spectrum is shown in Fig. 7, whereas the ratios 
between the ASCA and BeppoSAX spectra and the "com- 
posite" spectrum are shown in Fig. 6. Both these ratios axe 
consistent with 1 below 1 keV. Above 1 keV their behaviors 
differ: the ratio against the BeppoSAX spectrum decrease 
monotonically up to 10 koV; the ASCA ratio remains ba- 
sically constant around a value of about 0.6, with a local 
excess around the iron Kq iron line. 

In principle, applying this "composite" model to the 
ASCA and BeppoSAX spectra should provide the ultimate 
answer to the question on which source in the NGC 5643 
field is responsible for the flux and spectral variations ob- 
served across the X-ray history of NGC 5643. Unfortunately, 
the composite model has too many parameters, and would 
therefore largely overfit the ASCA and BeppoSAX spectra. 
We have therefore limited the scope of the exercise, keep- 
ing all the parameters of the composite model frozen to the 
XMM-Newton values, except the continua and iron line nor- 
malizations. This represents of course a very crude (and ad- 
mittedly arbitrary) assumption. The results of this exercise 
are summarized in Table 7. Its main conclusions can be out- 
lined as follows: 

• the ASCA soft X-ray faint state might be mainly due 
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Energy (keV) Energy (keV) 

Figure 6. Top: spectra {upper panels) and residuals in units of standard deviation {lower panels), when the ASCA {left) and Bep- 
poSAX {right) spectra are fit with a power-law model. Bottom: ratio between the ASCA and BeppoSAX spectra and the XMM-Newton 
"composite" model (details in text). 



Table 7. Continua and iron line normalizations when the "com- 
posite" model is applied to the ASCA and BeppoSAX spectra. 
The XMM-Newton best-fit values arc shown as well for compari- 
son. Legenda: Npi^- normalization of the locally absorbed power- 
law; Npi^^: normalization of the locally unabsorbcd power-law; 
^th,s - normalization of the softer thermal component; N^h f^: nor- 
malization of the harder thermal component. The superscripts " 
and refers to the NGC 5643 nucleus and sources X-1, respec- 
tively. Errors are at 1-a level for 1 interesting parameter. Values 
are expressed as xlO~^ of the corresponding units. Power-law 
norinalizat ions ar(* a1 1 kcV 



nuclear emission, leaving source X-1 as the dominant con- 
tributor to the hard X-ray flux 

It is worth noticing that the intensity of the iron line is 
consistent with being constant across all the observations. 

Given the crudoncss of the above assumption, these con- 
clusions must be regarded as no more than hints to the 
true physical picture. The ultimate answer on the nature 
of the X-ray variability in this system must await future 
re-observations with either Chandra or XMM-Newton. 



Normalization 


ASCA 


BeppoSAX 


XMM-Newton 


N" 

pi, a 


7.2 ± 1.9 


< 0.6 




N" 

pl,u 


ai+0.04 


< 0.5 


0.27 ±0.02 


pi, a 


< 0.4 


1.4 ±0.6 


1.97 ±0.05 


tn,s 


4±3 


< 18 




fjn 
"th,h 


0.7 ±0.2 




0.83 ± 0.05 


I" 

line 


1.8 ±0.3 


2.o±;:i 


'■■"^0.2 



to source NGC 5463 X-1 becoming at least five times fainter 
than during the XMM-Newton (or BeppoSAX) observation 
• the hard X-ray steep BeppoSAX spectrum might be 
mainly due to a decrease of the flux associated with the 



4 DISCUSSION 

4.1 The nature of the obscured AGN in 
NGC 5643 

A qualitative difference exists between objects obscured by 
matter with a column density smaller (Compton-thin) or 
larger (Compton-thick) then Nh — cr^^ — 1.5 x 10^'' cm~^ 
at energies lower than 10 keV, where most of the X-ray de- 
tectors flown so far were sensitive. Compton-thick AGN can 
be observed only along optical paths, which do not intercept 
the absorbing matter. This may make highly uncertain the 
determination of their intrinsic luminosity, which is depen- 
dent on the largely unknown distribution and covering frac- 
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Figure 7. XMM-Newton "composite" spectrum, combining 
NGC 5643 nucleus (Model #3 in Table 3), and X-1 source best- 
fits. Pale solid lines: thermal components; dashed lines: AGN 
continuum and Kq iron fluorescent line; dotted line: X-1 best-fit 
power-law; thick, upper solid line: total composite spectrum 

tion of the reflecting and/or scattering matter. Even more 
importantly, a Compton-thick absorber substantially sup- 
presses the incident radiation even at energies larger than 
the photoelectric cut-off. At 30 keV (100 keV), the fraction 
of transmitted radiation with respect to an absorber with 
Nh = 10^^ ''^ cm"^ is about 85, 30, 4 per cent (50, 20, 2 per 
cent) for Nh = lO^^ ^s cm^^, lO^^-^^ ^^^^-2^ ^q25.25 ^^^^-2^ 

respectively (Wilman & Fabian 1999). Knowing how many 
Compton-thick AGN exist may have an impact on the his- 
tory of accretion in the universe (Fabian 1999) 

Previous BeppoSAX observations suggested that 
Compton-thick objects represent a large fraction of the 
Seyfert 2 population (Maiolino et al. 1998), maybe as large 
as 30-50 per cent (Risaliti et al. 1999) . To identify Compton- 
thick AGN, the measurement of two critical observables 
was required: the shape of the power-law continuum above 
2 keV, and the EW of the iron line. Observed flat spec- 
tra (F < 1.0), and large EW (> a few hundreds eV) have 
been considered signatures of a reflection-dominated (i.e., 
Compton-thick) AGN. In principle, a detection by the Bep- 
poSAX PDS above 15 keV was crucial in determining the 
nature of the X-ray absorption (Matt et al. 1997b; Matt et 
al. 1999). However, the PDS sensitivity (a fraction of mCrab 
in a typical 50 ks observation) prevented robust detections 
in most Seyfert 2 galaxies. Based on the BeppoSAX re- 
sults, Maiohno et al. (1998) had classified the NGC 5643 nu- 
cleus as a "warm-scattered" Compton-thick AGN, obscured 
by Nh > 10^^ cm~^. This classification was mainly driven 
by the steep, unabsorbed spectrum observed by BeppoSAX, 
which is inconsistent with Compton-reflection dominance. 
It is interesting to observe that this classification was not 
supported by the measurement of the Kq, fluorescent iron 
line centroid, which was (and has been across the whole X- 
ray history of NGC 5643) inconsistent with He- or H-like 
iron fluorescent Ka transitions. XMM-Newton measured a 
Compton-thin absorber instead, with a column density in 
the range A'^ff6-10xl0^^ cm~^. The absorber either directly 
covers the nuclear emission, or its Compton-reflection. In 
the latter scenario, we might be observing reprocessing from 



the inner far side of the "torus", partly absorbed by the 
outer rim of its near side. This scenario would require a 
rather fine-tuned orientation. This, together with the lack 
of appropriate statistics, might explain why this kind of 
scenario has never been invoked so far in the context of 
reflection-dominated AGN. Incidentally, the lower column 
density helps to substantially reduce (but not to completely 
reconcile) the need for a dust composition geared toward 
large dust grains implied by X-ray obscuration with respect 
to the observed narrow Br^ fine component fluxes (Lutz et 
al. 2002). 

The lack of detection of transmitted nuclear emission in 
the Compton-reflection dominated scenario might be due to 
the line-of-sight toward the nucleus intercepting Compton- 
thick matter. This may indicate that the "torus" is thicker, 
the closer to its mid- plane one looks (Matt et al. 2000). Al- 
ternatively, the line of sight to the nucleus during the XMM- 
Newton observation could be probing a transient phase of 
low activity. The recent serendipitous discovery of transi- 
tions between transmission- to reflected-dominated spectra 
of Seyfert 2 galaxies (see Matt et al. 2003, and references 
therein) shows that the classiflcation of an obscured AGN as 
Compton-thick/-thin can be time-dependent. In at least two 
weU studied cases (NGC 2992; GiUi et al. 2000; NGC 6300, 
Guainazzi 2002), these transitions are due to changes of the 
overall AGN output by more than one order of magnitude, 
on timescales of the order of a few years. The reflection- 
dominated state in these cases is the "echo" of a previous 
state of AGN activity, as in the "swan song" state of the Nar- 
row Line Seyfert 1 galaxy NGC 4051 (Guainazzi et al. 1998; 
Uttley et al. 1999). Interestingly enough, earlier ASCA and 
BeppoSAX spectra of the NGC 5643 nucleus are markedly 
different in both flux and spectral shape from those observed 
by XMM-Newton. Taking into account the possible contam- 
ination of "source X-1" into the large ASCA and BeppoSAX 
aperture, the observed spectral variability is consistent with 
an historical dynamical range in the AGN power of one order 
of magnitude. Indeed, our interpretation of the steep spec- 
trum observed by BeppoSAX is that the nucleus was out- 
shone by NGC 5643 X-1, due to a phase of very low AGN 
brightness. This interesting possibility can be tested with 
future monitoring campaigns of this active nucleus, that we 
plan to pursue in the nearby future. Of course, we cannot 
rule out that the steep BeppoSAX spectrum is due to a re- 
vival of the AGN itself. This would make the X-ray history 
of NGC 5643 even more intriguing, though. 

The nature of the soft X-ray (i.e. < 2 keV) emission 
in the NGC 5643 nucleus is still uncertain. All possible 
models for the soft X-ray spectrum requires a contribution 
from emission lines. The quality of the data is, however, not 
good enough to establish the physical process responsible 
for them. A good fit is obtained with the composition of 
two thermal emission components from optically thin, coUi- 
sionally ionized plasma, with temperatures kT ~ 0.15 keV, 
and kT ~ 0.67 keV. A possible origin for this emission is gas 
heated by shocks or winds in regions of intense star forma- 
tion. The 0.5-4.5 keV luminosity associated with these ther- 
mal components (I/0.5-4.5 keV — 1-4 x 10*'' erg s~^) is con- 
sistent with the empirical relation discovered by David et al. 
(1992) with the Far InfraRed (FIR) luminosity in starburst 
galaxies; NGC 5643 Lfir = 1-2 x lO'^L© ~ L°o:f_^,s keV 
(Genzel et al. 1998; Moran et al. 2000). On the other hand. 
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the 2-10 keV luminosity ol the "ultra-hot" thermal com- 
ponent in Model#2 of Table 3 (~ 5 x 10^^ erg s"^) is 
only marginally consistent with the empirical relation be- 
tween 2-10 keV and FIR luminosity established by Ranalli 
ct al. (2003). Conversely, one can estimate the bolometric 
luminosity of the AGN in NGC 5643, L^„f ^ - 30Lx ^ 
7.5 X 10*^ erg s~^, if it is has a typical Spectral Energy 
Distribution as in, e.g., Elvis et al. (1994) and one assumes 
Model#3 in Table 3. A substantial contribution of AGN re- 
processing to the measured FIR luminosity would not be 
required. However, evidence for an intense nuclear starburst 
in this galaxy is still controversial. Alternatively, the EPIC 
soft X-ray spectrum can be fit with a pure photoionizcd 
plasma, whose signature in the optical band would be rep- 
resented by the bright one-sided ionization cone (Simpson et 
al. 1997). In this scenario, it is ruled out that a substantial 
contribution to the soft X-ray emission comes from scatter- 
ing of the nuclear continuum, as this would imply too steep a 
slope to be consistent with that inferred from the hard X-ray 
spectrum of this very same observation. However, it is pos- 
sible to fit the soft X-ray spectrum with a pure emission-line 
model. High-resolution spectroscopy observations of bright 
Seyfert 2 galaxies have indeed shown that the contribution 
of scattered continuum to the photoionizcd emission could 
be negligible (Kinkhabwala et al. 2002). In this scenario the 
line identification is not fully unambiguous, and contami- 
nation between transitions of different elements likely. At 
Oth-order, the interpretations fall into two broad classes. In 
the first the spectrum is dominated by He-like Ka transi- 
tions of Carbon, Oxygen, and Neon, with features from the 
last element being particulaxly abundant and including as 
well a Kff line, and a Radiative Recombination Continuum 
(RRC). Interestingly enough, the centroid of the feature en- 
compassing the Ovii triplet - which EPIC data cannot re- 
solve - is consistent with intercombination and resonant lines 
being the dominant contributors of the triplet, by contrast 
to NGC 1068. This suggests a large density (n > 10^^ cm"^) 
for the photoionizcd medium (Porquet & Dubau 2000). In a 
second possible class of interpretations, the line spectrum is 
dominated by transitions of highly ionized iron species from 
Fcxvil to Fexxiv. This interpretation would be consistent 
with the detection of a weak Fexxvi line in the hard 
part of the EPIC spectrum (cf. Fig. 4). It is interesting to 
observe that our line deconvolution of the EPIC spectra con- 
tains all the brightest lines discovered in the high-resolution 
spectrum of NGC 1068 (Kinkhabwala et al. 2002), except 
its Oviii and Nvii. Ovii/3 is missing as well, although in this 
case it might be blended to invisibility with the Nvii RRC. 



4.2 NGC 5643 X-1: an ultrabright ULX? 

The XMM-Newton image of the NGC 5643 field clearly 
shows a source, about 0.8' NE from the position of the 
optical nucleus. This source - labeled "NGC 5643 X-1" in 
this paper - is about a factor 50 per cent brighter than 
the nucleus itself, and had been previously observed by 
the ROSAT HRI, in a three times fainter state. In Fig. 8 
we show an overlap between the EPIC intensity contours 
and the simultaneous OM image. Although the relatively 
broad Point Spread Function of the XMM-Newton optics 
prevents us from pinpointing an unambiguous optical coun- 
terpart for NGC 5643 X-1, it is apparently located at the 
outskirts of the host galaxy optical extension. If the source 
is not a background object, its remarkable X-ray luminos- 
ity (~ 4.4 X 10"*" erg s^^) would convert it in the third 
brightest ULX ever observed in hard X-rays (Foschini et 
al. 2002; Humphrey et al. 2003; Swartz et al. 2003). Its 
spectral properties provide relatively little information on 
its nature. A fit with a simple power-law (F ~ 1.7) or 
bremsstrahlung (kT ~ 8 kcV) continuum with moderate 
absorption {Nh — 1 2x10'^^ cm^'^) is a good description of 
the observed spectrum. It is interesting to observe that a fit 
with a multitemperature disk blackbody is poor. 

In the quest for the nature of NGC 5643 X-1 the 
discovery of a possible counterpart at other wavelengths 
would be crucial. No known radio source can be asso- 
ciated with the position of NGC 5643 X-1. A 2MASS 
bright optical source with coordinates: 02000 = 14''32'"43.0''; 
^2000 ~ — 44°09'40", about a factor of 3 fainter than the 
optical nucleus, is only marginally consistent with typical 
uncertainties in XMM-Newton attitude reconstruction. The 
brightest OM UV source in the EPIC error box has co- 
ordinates a2()()() = 14''32'"42.9''; ^2000 = -44°09'34". The 
identification of an unambiguous counterpart needs to await 
a more precise determination of the X-ray centroid, which 
could be achieved with a Chandra observation of the galaxy 
field. On the other hand, the study of the X-ray history of 
the NGC 5643 field suggest that NGC 5643 X-1 underwent 
flux variations with a dynamical range of at least 5, thus 
ruling out an association with a supernova remnant. Again, 
future monitoring campaigns of this field may offer impor- 
tant clues on its variability pattern, and henceforth on its 
nature. 



The quest for the physical origin of the soft X-ray emis- 
sion in this object requires instrumentation of enough power 
to be able to perform high-resolution spectroscopy in the 
soft X-ray regime. Chandra high-resolution imaging could 
be important as well in determining whether the soft X-ray 
nuclear spectrum is dominated by diffuse emission, associ- 
ated to the 15" ionization cone. AGN photoionizcd plasma 
structures extending on scales of hundreds parsecs have been 
observed in nearby Seyfert 2s (Sako et al. 2000; Bianchi et 
al. 2003; Iwasawa et al. 2003). Evidence from the XMM- 
Newton image is still inconclusive in this respect. 
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Figure 8. Left panel EPIC pn iso-intcnsity contours superposed to the innermost science window of tlie OM/UVWl exposure. The 
contours represent 9 logarithmically scaled steps in the pixel count range between 5 and 169. Right panel: zoom of the innermost 1'50". 
The pn iso-intensity contours are superposed to the K-band 2MASS image 
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